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O ■ Abstract 

We present a quenched lattice calculation for the lowest lying 66g-hybrid states in the 
framework of NRQCD using the leading order Hamiltonian up to 0(mv 2 ). We demon- 
strate the existence of a nearly degenerate rotational band of states with an excitation 
| \ energy approximately 1.6 GeV above the T ground state. This lies around the BB* r 

■ threshold but well above the £?_B-threshold. Therefore a heavy hybrid signal may well 
(■-<■) \ be detected if the centre-of-mass energy in B-factories is raised a few hundred MeV to 

coincide with other resonances above the 4S state. Our prediction is consistent with most 
phcnomcnological models and lattice calculations carried out in the static limit. 

'. 1 Introduction 

Hybrid mesons are of intense interest both theoretically and experimentally because of the 
opportunity they provide for investigating nonperturbatively the gluonic degrees of freedom in 
QCD. In contrast to the standard qq- mesons in which the quarks form a colour singlet, hybrids 

■ contain quarks in a colour octet state. There exist hybrid states that have quantum numbers 
J> \ not available to pure qq-st&tes and as a result do not mix with them. These exotic states are 

of particular interest. Interest has been further heightened by the recently reported discovery 
at Brookhaven of a 1 h state at (1370 ±16 MeV [||. However, the experimental study 
of such light hybrids is made difficult by the density of levels in the 1-2 GeV range and by 
strong mixing effects. Such difficulties are minimised for heavy hybrids and this should result 
in a clearer signal at the appropriate energies in future B-factories. 

The study of hybrid states has been approached in a number of ways, such as flux tube 
models bag models Q, sum rules § and the constituent gluon model ||. Light hybrids 
have been studied in the framework of lattice QCD || |7| . Heavy hybrids have been approached 
through the static quark limit § by identifying a (^-potential appropriate to each gluonic 
excitation ||. 

In this paper we study hybrid excitations of the 66-system using NRQCD |Tot| . This allows 
us to go beyond the static limit. Our calculation is preliminary in that we use only the lowest 
order heavy quark Hamiltonian and employ the quenched approximation for the gluon field. 
Previous work confirms that this approximation gives a reasonable account of low lying spin 
averaged energy levels of heavy quark systems [11, 12]. In this sense we believe we have 
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achieved an acceptable computation of the lowest magnetic hybrid excitations. We have not 
yet computed the hybrid levels with opposite parity. Within bag models these electric hybrids 
are believed to lie higher in energy than the corresponding magnetic states Q. 

In Section 2 we introduce the non-relativistic evolution equation for the quark propagator. 
In Section 3 we discuss continuum operators that connect the hybrid states to the vacuum and 
in Section 4 we set out the lattice versions of these operators that were used in our simulation. 
The results are presented in Section 5. 



2 NRQCD and Heavy Quark Propagators 

The NRQCD approach to the computation of the heavy quark propagator G(x,y), has been 



explained previously [11, 12 1. For completeness we record here the Euclidean time evolution 
equation, namely 

G(x,t + l; y ) = (l-^)"[/t(x)(l-^) n G(x,t;y) , t > t y . (1) 

The initial condition has the form 

G(x,t = V,2/) = S(x,y) . (2) 

where S(x,y) is the source term on the first timeslice, (t = t y ), appropriate to the channel 
under study. In this paper the Hamiltonian is 

A 2 

Ho = -^~ , (3) 
where mj, is the bare quark mass and A 2 is the standard spatial covariant Laplacian defined 



in [10]. In the calculation all the link variables are tadpole improved according to the replace- 
ment Ufj,{x) — * U fl (x)/uo where no is obtained from the plaquette Uu- uo = (0|^TrC/n I0) 1 / 4 . 
Other suggested improvements [13], make no significant difference to the present calculation 
which is concerned only with spin averaged quantities. This is consistent with our neglect 
of higher corrections to the Hamiltonian that incorporate quark spin-gluon couplings and 
C(mf 4 ) terms. 



3 Hybrid States and Operators 



We will use the standard nomenclature for hybrid states [14] in which charge conjugation and 
parity satisfy 

C = (-l) l + s+1 (4) 

P - { { ~ 1)l+J TE (5) 
\ (-1)^' +1 TM ' W 

where I and s are the qq orbital angular momentum and spin, and j is the gluon angular 
momentum. The historical notations TE and TM refer to the magnetic and electric hybrid 
states, respectively. 

To extract masses we calculate two-point functions of operators with the appropriate 
quantum numbers. For hybrid states these operators must include a gluon field factor in 



2 



order to ensure the presence of a gluon excitation. This is done by introducing the gauge field 
operators 

Bi = -j^ijkFjk , (6) 
Ei = F it , 

where is the gluon field tensor. We can alternatively replace E by V x B to access the 
same quantum numbers on a single timeslice. 

From Bi we can construct the j = 1 magnetic hybrid operators. The spin singlet states 
are coupled to the vacuum by the operators defined in Table [l]. 



state 


I 


J PC 


X^BiiP 





\— 


xHBi,DiW 


1 


0++ 


X^eijk{Bj,D k }il) 


1 


1++ 


xH{Bi,Dj} + {Bj, A} - iSijiB^DjW 


1 


2++ 



Table 1: Continuum operators for spin-singlet hybrid states. The non-relativistic 2 spinors, 
tp and x j denote the quark field and anti-quark field, respectively. 

It has been pointed out by Griffiths et al. that there is a common gg-potential in- 
teraction for all the j ' = 1 states in the static quark limit. The presence of excited glue, 
however, results in a relatively shallow structure for the radial dependence of the potential 
function. Consequently we can anticipate that the orbital motion of finite mass quarks in the 
presence of this potential will give rise to a nearly degenerate rotational band of states. In 
particular the I = and I = 1 states will be approximately degenerate. If, in addition, the 
orbital dynamics of the quarks is controlled by a spin-independent force, as is the case in the 
approximation used in this paper, the triplet quark states built on the above singlet states 
will also lie in the same degenerate band. The triplet states shown in Table |2| are built on 
the 1 state and contain the 1 ^ which has exotic quantum numbers. Triplet states built 



state 


I 




X'o-iBilp 





o-+ 


eijkX^jBkip 





1-+ 


xH a iBj + OjBi - ^Sij(7iBi)ip 





2-+ 



Table 2: Continuum operators for spin-triplet states corresponding to the 1 

on the 1 ++ are shown in Table ||. Those contain also the exotics 0" 1 and 2^ . There are 
further triplet states including a 3" 1 state. All these states can be expected to be essentially 
degenerate in the case of our simple spin-independent Hamiltonian. 
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state 


I 


J P C 




1 


o+- 


X^jd^i, Dj} — {Bj, Di})ip 


1 


1+- 


X ] {vi(jki{B k ,Di} + aje ik i{B k ,Di} - ±5ij(Jie jkl {B k , Di})ip 


1 


2+- 



Table 3: Continuum operators for spin-triplet states corresponding to the 1 ++ . 

4 Lattice Operators 

Because our evolution equation involves no spin-corrections to the Hamiltonian we focus on 
the lattice versions of the operators in Table |]. In constructing them we replace covariant 
derivatives, Di, with covariant lattice derivatives, A", in the extended form used in [12, 16]. 
For the magnetic gluon operators which involve only spatial derivatives this presents no 
problem since these operators can be formulated in terms of variables on a single timeslice. 
We define the colour magnetic field on the lattice as 

£r = ^ fc [A™,A£] , (7) 

with the definition 



A^(x) = Uli){x + m) - L n _${x - ni) 

Lfip(x) = Ui{x)Ui{x + i)...Ui{x + {n-l)i)i){x + ni) (8) 

As in an earlier study |l6| , these operators result in commutators of extended link variables 
for the hybrid state 

Bt = \e ijk {[L],Ll]-\L- p Ul k ]} , (9) 

For the free field case, where U^{x) = 1, these operators vanish as expected. Finally we 
have the extended hybrid operator 

fl?(x)= 6 y jkX t (x)([Ly,L3f]-[L» J ,L2 fc ])^(x) . (10) 

We have achieved a significant improvement for the signal by employing the fuzzing al- 
gorithm for link variables suggested in We use a central link weight of c = 2.5 with six 
fuzzing iterations. We now use those fuzzed link variables to construct the extended links, 
which we then use in the meson operators. Phenomenological models indicate that hybrid 
states are more extended than the standard (/(/-states. This suggests that operators with large 
spatial separation will have the best overlaps. In compromising between spatial size of the 
operators and the efficiency of the code we choose n=4 and 5. In some cases the operators 
were further improved by using, in addition, Jacobi-smearing for the quark fields flis} |. 

The meson correlator is written as a Monte Carlo average over all configurations 

C nm (x, y) = (tr \G\x, y)G nm (x, y)] ) , (11) 
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O h irrep. 


xHx) 


lowest contiuum J 


M 


1 


o-+ 


Ti 


eijk AjAfe 


1— 


Ax 


eijkAiAjAk 


0++ 


T x 


Cijk £fcim{Aj, A;A TO } 


1++ 


T 2 


Sijk£klm{Aj, A;A m } 


2++ 


E 


S a jk£klm{Aj, A;A m } 


2++ 



Table 4: Equivalent lattice operators of Table |l]. The first column denotes the irreducible 
representation of the octaheadral group. We define Sijk = \tijk\ and S a jk projects out the two 
linearly independent traceless-symmetric combinations corresponding to the representation 
E. 



p 


6.0 


L 3 x T 


16 3 x 48 


a- 1 in GeV (from IP - IS) 


2.44(4) 


arrih, n 


1.71, 2 


u 


0.878 


spatial starts 


8 


temporal starts 


5 


total number of measurements 


20, 000 



Table 5: The parameters in our simulation. To increase the statistics, we chose several starting 
points per configuration. 

where tr denotes contraction over all internal degrees of freedom and Q nm is the smeared 
propagator defined by 

G nm (x,y) = £ O n (x,z 1 )G(z 1 ,z 2 )0^(z 2 ,y) . (12) 

21,22 

Here (n, m) stands for the radii at (sink, source) and the O n are operators as defined in Table 
[|. For the extended propagator we solve equation [l] with S*(x, y) = O m ^(x, y) and multiply 
with O n at the sink. We fix the origin at some (arbitrary) lattice point, y, and sum over all 
spatial x so as to project out the zero momentum mode. 

5 Simulation and Results 

The parameters of our simulation are shown in Table || The quenched gauge field configura- 
tions were all generated at the EPCC in Edinburgh. The propagators were calculated at the 
HPCF in Cambridge. 

We fit the correlators to the multi-exponential form 

Cnt) = Y,<Te- M ^ ■ (13) 
1=1 
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Operator 


State 


Fit r 

a^ 


3SultS 

aEi — aE 


Energies 

Eq — E~£ 


\ [GeV] 
Ei — E 


1 


o-+ 


0.4487(13) 




0. 






1" 


1.114(40) 


0.628(66) 


1.62(10) 


1.53(16) 


A* • Bi 

{A { i,B 2} } 
{Ai,Bi}-{A 2j 5 2 } 


0++ 

1++ 

2++(T 2 ) 
2 ++ (E) 


1.175(22) 
1.134(19) 
1.161(19) 
1.126(13) 


0.609(73) 
0.61(10) 
0.600(61) 
0.520(58) 


1.775(61) 
1.674(54) 
1.740(55) 
1.655(42) 


1.49(18) 
1.49(25) 
1.47(15) 
1.27(14) 



Table 6: Hybrid spectrum from 499 configurations at /3 = 6.0, am^ = 1.71, no = 0.878, 
a lP-iS = 2.44(4) GeV. Eq and E\ denote the energies of ground and first excited states in 
the relevant channel. The results are obtained from a two-exponential fit as shown in column 
3 and 4. In column 5 and 6 we have converted all results into physical units, using the inverse 
lattice spacing which is determined from the IP — IS splitting. The ground state in each 
channel is given relative to the T(15). 



Here a denotes a meson state with certain quantum numbers, (n, m) the different radii at 
the (sink, source) and t the Euclidean time. As the hybrid excitations are very high above 
the ground state the signal-to-noise ratio deteriorates very quickly and vanishes at t £ 15. 
Therefore we choose i max = 12 and vary t m \ n in search of a plateau. An example of such a 
plot for the lowest lying bb- hybrid is shown in Figure [l|. A collection of our results can be 
found in Table |6| and Figure |2| 

We quote the masses from plateaus where x 2 /dof < 1. The error estimate is calculated 
from the inverse of the Hessian matrix which is determined during a Levenberg-Marquardt 
X 2 minimisation [19]. 



6 Conclusions 

In this paper we reported on a heavy hybrid signal obtained within the framework of NRQCD. 
This field theoretical calculation goes beyond the static approximation and so has the virtue of 
incorporating the dynamics of the heavy quarks. In our investigation we calculated the masses 
of hybrids where the (/(/-pair is in a colour octet coupled to the colour magnetic field. Some of 
the results have already been reported elsewhere |2C|j . Work on the same configurations using 
different operators has been discussed in |Hf| . We obtained a signal for several correlation 
functions where the magnetic field couples to spin-singlet states of different orbital angular 
momenta, as listed in Table |j. Static potential models predict that the above states should 
lie in a nearly degenerate rotational band. From Table || and Figure |2| it can be seen that our 
results confirm this picture to within two standard deviations. For our Hamiltonian all the 
spin-triplet states are degenerate in energy with the spin-singlets and we are only sensitive to 
spin-averaged quantities. Therefore, the exotic 1 *~ state is degenerate in this approximation 
with the 1 and likewise both the 0" 1 and 2" 1 are degenerate with the 1 ++ . 

Phenomenologically, there are two important thresholds for the production and decay of 
hybrid ^-states. These are the 5I?-threshold at 10.56 GeV and the B Bj-threshold at 11.01 
GeV. For the latter we assume that B*j (5732) is indeed a P-wave as suggested in the Particle 
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1— exponential fit 
o 2— exponential fit 
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Figure 1: Fit to hybrid signal. Consistent fit results are shown for the ground state for both 
a 1-exponential fit and a 2-exponential fit. 
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Figure 2: Results. The 1 is 1.62(10) GeV above the IS. The corresponding rotational band 
appears to be degenerate around the i?i?}-threshold. Each state is the spin-averaged and 
some of them contain also the exotics 1~ + , + ~, 2 + ~. See main text. 
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Data Book |2^| . Sometimes this state is also refered to as B** . Below the (S+P)-threshold 
hybrid states are thought to be stable [Q. As can be seen from Figure |2| our results suggest 
that hybrid states lie close to the (S+P)-threshold and 4-5 standard deviations above the 
(S+S)-threshold. 

In our calculation there are a number of sources of systematic errors. First, we have 
retained only terms 0(mv 2 ) in the Hamiltonian and neglected those 0(mv 4 ) and higher. 
Numerically this may not be a bad approximation since we expect the quarks in the hybrid 
states to be even more non-relativistic than in the T itself. This is consistent with the 
very shallow qg-potential predicted in the static limit of the hybrid state H, ^4], and the 
implied near degeneracy of the resulting rotational band - a result confirmed in our calculation. 
However a study of singlet-triplet splitting requires that higher order terms, such as a ■ B, be 
taken into account. The above argument suggests that this splitting is small. Furthermore, 
there may be finite size effects because the hybrids in the shallow potential are expected to 
be larger than the T. From studies in the static limit one expects the interquark seperation 



in bbg to be of the order of 0.5 fm [21, 25|. The actual extent of a hybrid state may still be 
bigger than this ||. We note that for the parameters used in this paper the lattice has a 
spatial extent of approximately 1.3 fm. 

Finally, and perhaps most importantly, we must take into account uncertainties in the 
value of a -1 , the inverse lattice spacing. These uncertainties are intrinsic to any calculation 
based on a quenched gluon approximation, since not all mass ratios can be simultaneously 
correct. The result therefore depends on which observables are used to fit a -1 . In the numbers 
quoted above we used the IP — IS mass difference of a standard NRQCD calculation of the 
T-system. This is a consistent approach and yields a -1 = 2.44(4) GeV and our number of 
11.08(10) GeV for the lowest hybrid. Perantonis and Michael quote 10.81(25) GeV from a 
Schrodinger Equation using the static potential from the lattice, where a" 1 is determined from 
the string tension in the quenched approximation [Q. Their error also includes an estimate 
for quenching effects. If we use their value of a" 1 =2.04(2) GeV, we obtain the mass of the 
1 h to be 10.82(8) GeV, which is consistent with their result. These uncertainties should be 
resolved in a simulation with dynamical light quarks where the coupling runs appropriately 
and the inverse lattice spacing is expected to be the same when determined at different scales. 
Therefore the issue whether the hybrid states lie above or below the i?i?*-threshold, must 
be addressed using unquenched gauge field configurations. Subject to this uncertainty in the 
value of the inverse lattice spacing, our results are consistent with the predictions from most 
phenomenological models [p6f |. 

Acknowledgements 



We would like to thank other members of the UKQCD collaboration in particular S. Collins 
for useful discussions. T.M. is supported by a grant from EPSRC (Ref. No. 94007885). 
H.P.S. is supported by the Leverhulme Trust. Our calculations are performed on the Hitachi 
SR2001 located at the University of Cambridge High Performance Computing Facility and 
the CRAY-T3D at the Edinburgh Parallel Computing Centre at Edinburgh University. 

References 

[1] D.R. Thompson et al. Phys. Rev. Lett. 79, (1997), 1630-1633. 



8 



N. Isgur and J. Paton, Phys. Rev. D31, (1985), 2910-2929. 

P. Hasenfratz et al, Phys. Lett. B95, (1981), 299-305. 
T. Barnes et al, Nucl. Phys. B224, (1983), 241-264. 

J. Govaerts et al, Nucl. Phys. B284, (1987), 674-689. 

Yu.S. Kalashnikova and Yu.B. Yufryakov, Phys.Atom.Nucl.60, (1997), 307-313 

P. Lacock et al, UKQCD Collaboration, Phys. Rev. D54, (1996), 6997-7009. 

C. Bernard et. al, MILC Collaboration, Nucl. Phys. B (Proc. Suppl.) 53, (1997), 228- 
230. 

C. Bernard et. al, MILC Collaboration, |hep-lat/9707008| . 

S. Perantonis and C. Michael, Nucl. Phys. B347, (1990), 854-868. 

I.J. Ford et al, Phys. Lett. B208, (1988), 286-290. 

CP. Lepage et al, Phys. Rev. D46, (1992), 4052-4067. 

C.T.H. Davies et al Phys. Rev. D50, (1994), 6963-6977. 

T. Manke et al, UKQCD Collaboration, Phys. Lett. B408, (1997), 308-314. 

H.D. Trottier, Phys. Rev. D55, (1997), 6844-6851. 

A. Le Yaouanc et al, Z. Phys. C28, (1985), 309-315. 

L.A. Griffiths et al, Phys. Lett. B129, (1983), 351-356. 

S. Catterall et al, UKQCD Collaboration, Phys. Lett. B300, (1993), 393-399. 

P. Lacock et al, UKQCD Collaboration, Phys. Rev. D51, (1995), 6403-6410. 

C.R. Allton et al, UKQCD Collaboration, Phys. Rev. D47, (1993), 5128-5137. 

W.H. Press et al, "Numerical recipes in Fortran", (2nd Edition), Cambridge University 
Press (1992). 

T. Manke et al, UKQCD Collaboration, to appear in Nucl. Phys. B (Proc. Suppl), 



proceedings of Lattice '97, Edinburgh, hep-lat/9709001 



S. Collins et al, UKQCD Collaboration, to appear in Nucl. Phys. B (Proc. Suppl), 
proceedings of Lattice '97, Edinburgh, |hep-lat/9710058 . 



R.M. Barnett et al, Phys. Rev. D54, (1996), 1. 

F. E. Close and P. R. Page, Nucl. Phys. B443, (1995), 233-254. 

K.J. Juge et al, to appear in Nucl. Phys. B (Proc. Suppl), proceedings of Lattice '97, 



Edinburgh, |hep-lat/9709132 



K.J. Juge et al, to appear in Nucl. Phys. B (Proc. Suppl), proceedings of Lattice '97, 
Edinburgh, [hep-lat/9709131. 



T. Barnes et al, Phys. Rev. D52, (1995), 5242-5256. 



9 



